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INTRODUCTION 


The results presented in Part II! of this final technical report 
are based on the conception of a Primary Tectonic and Seismic Motivating Force 
(PTSMF) and a procedure to test its verification by numerical calculation, both 
of which were conceived and developed prior to the initiation of the present 
contract. Consequently, this report pertains primarily to the numerical imple- 
mentation of this procedure, with the principal objective of numerically verify~ 
ing the PTSMF by calculating the stress and deformation fields which it produces 
in a hierarchy of mechanical models of the actual earth portraying its mechanic-~ 
al non- unifromities with increasing specificity. However, to give proper per- 
spective to this task, we present in Part | the historical development of the 
work which led to the conception of the PTSMF and in Part I| an outline of the 
theoretical considerations which help to show that the PTSMF was discovered in 
a conceptual framework which, in fact, rests on the bedrock of physical theory, 
and which may serve to demonstrate the primacy of force in relation to the priti- 
cioles of classical mechanics and physical theory as a whole. 

Although the numerical results presented here substantially verify this 
force as the principal motivator of tectonic and seismic phenomena in the earth, 
its scientific and practical ramifications remain to be fully explored in detail, 
which, unfortunately, we have been unable to do under the present contract due 
to insufficient funding. When the above task is completed in subsequent inves- 
tigations, we then plan to present in detail the conceptual, theoretical and an- 
alytical background of PTSMF that were developed prior to the initiation of the 
present contract, and which led to the numerical analysis and calculations con- 


ducted with its support. 


HISTORICAL BACKGROUND AND DEVELOPMENT 


For the purpose of rendering a reasonably self-contained final technical ~e- 


port and a coherent presentation of the unfolding of ideas that led to the numer- 
ical results presented, we shall recapitulate and reference some earlier ideas 
and their historical development, as partially reported in publications and lec~ 


tures, in this part of the report. 


A. A Phenomenological Mechanical Description of the Seismic Process and the 
Conception of a Phenomenological] Mechanical Model of the Earth 


The point of origin which led originally to the recognition that relax- 
ation phenomena in the earth are an essential and, in some cases, a dominant as- 
pect of important geomechanical processes, was the search for corceivable and 
physically edmissible processes which may account for the production and recur- 
rence of earthquakes in seismically active regions , where extra terrestial ag- 
ents and catastrophic forces suddenly produced in the earth are excluded as cas- 
ual agents. All cases considered, namely shallow, intermediate and deep focus 
earthquakes, led necessarily to the conclusion that relaxation phenomena in the 
materials of the earth are essential to account for these various kinds of earth- 
quakes. These considerations, in turn, led to the very significant and apparent-~ 
ly new conclusion, that the inelasticity of the earth is a fundamental aspect of 
its seismicity, an aspect which in fact was previously de-emphasized and therefore 
overlooked in seismological models and in their mathematical applications to 
theoretical seismology. 

The important inference cited above concerning the fundamental role of relax- 


ation phenomena in the mechanical evolution of earthquakes and by which it was 


concluded that the inelasvicity of the earth is an essential feature of its seis- 
micity, was originally inferred from a phenomenological description of the se.s- 
mic process which evidently incorporates indisputable phenomenological features 
of this peck? By seismic process or seismic mechanism we mean here a meciian- 
ism which (1) progressively accumulates strain energy in the crust, (2) redis- 
tribution of the accumulated strain energy in the crust by relaxation processes, 
(3) the consequent localization and elevation of stresses at the unrelaxed re- 
gions within the crust, (4) the dynamical release of these localized and elevated 
stresses when the crustal material faiis to support them in static equilibrium 
which phenomenologically defines the onset of an earthquake, and dynanically 
loads the sub-crustal material, (5) recovery of static equilibrium, phenomenolog- 
ically defined by the termination of the earthquake, including aftershocks that 
follow, and the re-emergence of queiscent state, (6) relaxation of stresses in the 
sub-crustal layer with consequent reaccumulation, redistribution, localization and 
elevation of stresses in the crust, and (7) the recurrence of earthquakes at the 
sites of elevated stresses, produced at newly localized regions in materials which 
can no longer support them in static equilibrium. This essentially completes the 
phenomenological description of the seismic process noted above which portrays 
the cyclicity of seismic phenomena in seismically active regions and in which 
certain critical levels of ambient stresses are produced and maintained by an 
appropriate seismic driving force. 

If, in addition, we postulate that the driving force depends strictly on 
location and is essentially time independent, then for reasons given above, the 


phenomenological description of the seismic process leads to certain definitive 


as eee ee 


conclusions concerning the mechanical phenomenological features of the earth which 
follow logically from the phenomenological description of the seismic process pre- 
sented and to which we must therfore assign the same credibility as we do the phe- 
nomenological facts on which the description is based. With this condition on the 
driving force, it follows directly from the above phenomenological description of 
the seismic process that the earth must necessarily be (1) non-uniform in mechan- 
ical constitution, (2) inelastic, and (3) at least partially Maxwellian inelastic. 
lt is noted that these basic inferences drawn logically from the phenomenological 


description are essentially phenomenological in character, as they should be, since 


they were inferred from a phenomenological description. In particular, the class 
of Maxwellian materials of the earth referred to in item (3) are phenomenological i 
rather than constitutively defined by the description of their behavior under phe- 
nomenologically prescribed mechanical conditions. Accordingly, the class of Max- 
welliar materials so defined, which includes as its simplest member the Maxwell 
body, in facr includes materials which constitutively have complicated time rate 
dependent and non-linear properties. Moreover, it is important to emphasize here 
that because these essential phenomenological ly mechanical features of the earth, 
which define a phenomenological mechanical model of the crust and mantle, are log- 
ically implied by phenomenological features of the seismic process which appear in- 
disputable, it is consequently expected that all subsequent findings in Geomech- 
anics should necessarily conform with this basic phenomenological mechanical mode] 


of the crust and the mantle. This expectation has, in fact, been consistantly sub- 


Stantiated by subsequent developments in Geomechanics up to the present time. 


id 


Ramifications and Specific Consequences of the Phenomenological 
Mechanical Model 


From our phenomenological description of the seismic process and the phe- 
nomenological mechanical model of the crust and mantle based upon it, there orig- 
inated a comprehensive study in geomechanics which has been sustained since their 
inception in 1950. 

(2) 

Our first attempt in 1952 was to explore analytically, in terms of compar- 
atively simple mathematical constitutive models, the implications of incorporat- 
ing Maxwellian visco-elastic properties in seismological models of the earth in 
order to assess their consequences on seismic wave propogation and observational 
seismology. These early investigations included the analytical study of the pro- 
pogation of Love waves in elastic crustal layers supported by a half space con- 
sisting of a Maxwell material. These led to stringent conditions on the admissi- 
bility of Love waves which are essentially determined by the bourdary conditions 
at the interface joining the elastic crustal layer with the underlying Maxwell 
half space. The stringency of these admissibility conditions implied a cut off 
of admissible periods and thereby indicated that the Maxwellian properties of the 
upper mantle, inferred from the phenomenological description of the seismic pro- 
cess, would produce, in seismic Signatures of earthquakes, features by which its 
specific inelastic mechanical properties can be identified. To our knowledge, 
this indeed was the first step taken, with purpose, toward introducing inelastic 
properties in mathematical seismological models of the earth and, in particular, 
Maxwellian inelastic features which had been inferred as essential for sustaining 


the seismicity of the earth. 


This investigation concerning Love waves was then extended to the study of 


Rayleigh surface waves in a half space consisting of a Maxwell material with the 


(3,4) 


purpose of establishing the conditions for the admissibility of such Rayleigh waves. 


As in the case of the surface Love waves, we also found here that the conditions 
for the admissibility of Rayleigh waves propogating over a Naxwell half space are 
more stringent than are the admissibility conditions obtained for Rayleigh waves 
propogating over an elastic half space. This resu’t further reinforced the con- 
clusion that the Maxwellian constitution cf the earth, and of the upper mantle in 
particular, should lead to discernable features in seismic signatures of earth- 
quakes, and conversely that these discernable features may be useful for indentify- 
ing the Maxwellian inelastic properties of the earth in terms of appropriate seis- 
mic models of the earth which mathematically represent its Maxwellian constitutive 
features. 

These two preliminary results gave good reason to conclude that the then pre- 
vailing tradition in tneoretical seismology which clearly overlooked and consequent- 
ly strongly de-emphasized the fundamental role of the inelasticity of the earth in 
producing and in sustaining its seismicity and its tectonic phenomena must be rad- 
ically changed. Indeed, had this fundamental fact been previously recognized, 
clearly every effort would have been made to appropriately incorporate, with par- 
ticular emphasis, in the mathematical seismological models, the inelasticity as a 
fundamental feature of the earth. With this realization, a procedure was conceiv- 
ed and then chee aon the elastic moduli appearing in various traditional 
mathematical seismological models of the earth are replaced by complex moduli. In 


this regard, we then emphasized that the pure imaginary parts of these complex 


moduli that represent stresses proportional to corresponding elastic Stresses but, 


which by their construction, lag them by 90° phase, are consequently non-conserv- 

ative, furnish a first approximation for representing the disipative effects pro- 

duced by very small deviations from perfect elasticity, without however represent~ 
ing in constitutive terms the nature of the inelasticity of the materials in which 
the dissipation is produced. 

Although it was understood at that time that this procedure provides a frame- 
work for relating wave attenuation to dissipation due to inelastic properties of 
the materials of the earth, it was also stressed that this model does not and can 
not give any information concerning the specific inelastic constitutive properties 
of the earth. This should be evident since one has to formally portray in the math- 
ematical models the essential constitutive features of the earth in order to deter- 
mine them numerically from seismic data. 

The above idea and procedure were communicated to D. L. Anderson in 1962, 
both in lecture and privately, and subsequently used by him and Archambeau in de- 
veloping a mathematical seismological model of the earth in terms of which they 
calculated from seismic data a distribution of the so-called quality index 'Q' 
in the earth. This was done apparently without realizing the importance of taking 
into consideration the nature of its inelasticity of the earth, and of its proper 
representation in constitutive terms as required for analytical and numerical in- 
vestigation of its observable and measureable consequences. Therefore, they er- 
ronously ccacluded that this quality index 'Q' is frequency independent, an error 
which can be directly attributed to the fact that their model did not include an 


equally important basic feature deduced from the phenomenological description of 


the seismic process, namely, its Maxwellian inelastic property as discussed earlier. 


This serious omission also led Anderson to the denial of the possibility of the 
(6) 


convection current in the upper martle. 


In this regard, it is instructive to note that recently a comprehensive 
paper by D. Dd. aul on attenuation of seismic waves emphatically reports that 
a critical examination of iia ws ea data in the context of a recent enalytical 
development of Backus and Kites ele inescapably to the conclusion that the 
thesis advanced by Anderson and ihcshal ce os the frequency independence of the 
'Q' of the earth is in fact untenable; a conclusion which follows directly from 
the phenomenological Maxwellian model of the earth originally inferred from the 
phenomenological description of the seismic process, and which had been reported 
by the senior author with considerable emphasis since :950. This is particularly 
noteworthy since D. D. Jackson had collaborated with D. L. eee the publi- 
cation of an earlier paper pertaining to the representation of the inelastic 
features of the earth in terms of the quality index 'Q', and he is consequently 
well versed with the basis of the inferences drawn by Anderson and Archambeau on 
the frequency independence of 'Q'. Moreove., the expression for the internal fric- 
tion presented by Jackson on page 31 of his paper, for the accomodation of atten- 
uation data and which clearly displays its dependence on frequency, does in fact 
correspond to a Maxwellian material. 

This development significantly substantiates the thesis that all subsequent 
developments in qeomechanics which are compatible with reliable measurements must 


necessarily comply with the three essential mechanical phenomenological features 


of the phenomenologicai mechanical model of the earth which prove to be central in 


the dynamical processes that attend the evolution of earthquakes ascribed to the 


primary seismic motivating force as conceived later in 1965. Even though we were 


in fact responsible for the idea to introduce complex moduli as a first step for 
extending the mathematical theories of elastic seismology to corresponding theor- 
ies of inelastic seismology, we were nevertheless well aware that this approach 
was essentially limited, particularly since such mathematical theories will not 
circumscribe certain essential features of the phenomenological mechanical niodel. 
This led to the second stage of our work which concerns the development of ade- 
quate and admissible mathematical theories of inelastic seismology which, for rea- 
sons already given above, must necessarily agree with the essential features of our 
fundamental phenomenological mechanical model of the crust and mantle. The devel- 
opment of such admissible theories of inelastic seismology accordingly demands the 
reconciliation of constitutive Maxwellian models of the upper mantle with seismic 
data. 

This second stage was initiated by modeling the crust and upper mantle as a 
layered configuaration, each layer consisting of the simplest Maxwellian material, 
namely a Maxwell body constitutively uevings As a crucial test of this model, we 
inquired analytically concerning conditions for its reconciliation with Love wave 
attenuation data. This reconciliation, however, demands that the viscosity co- 
efficient of the layers representing the upper mantle in this model be many orders 
of magnitude smaller than the accepted value calculated by wBeeeTT using the Fen- 


noscandian recoil data following the removal of the iceload. Clearly, such recon- 


wpe 


ciliation would be untenable if Haskell's calculated value is indeed representative 


of the actual viscosity of the upper mantle, and this, of course, would seriously 


challenge the credibility of our basic phenomenological mechanical model. 


However, since the phenomenological mechanical model was logically inferred 


from the phenomenologica! description of seismic processes, which includes virtu- 
ally indisputable facts, our confidence in the validity of this phenomenological 
mechanical model was naturally sustained and this gave compelling reason to crit- 
ically examine the theoretical basis for Haskell's calculations. This examination 
led to the cruciai question; whether the Fennoscandia recoi] data which Haskell 
used is reconcilable with a three layered half space of Maxwellian material which 
includes a thin layer of very low viscosity situated under the crust. This ques- 
(14,15) 
tion was formulated and resolved affirmatively in mathematical terms. In addition, 
Haskell's theo-y is recovered in this formulation as a limiting case when the two 
upper layers are removed. Therefore, this basic phenomenological mechanical model 
and this particular constitutive model based upon it readily accomodate both the 
se smic phenomena which are produced in a very short time scale as well as the re- 
coil data associated with the Fennoscandia uplift that occured over a very long 
geological time period. This joint accomodation demands the existence of a thin 
and very low viscosity Maxwellian layer underlying the crust. It should be evi- 
dent that this has important ramifications regarding continental drift, convect- 
ion currents, plate tectonics, etc. Moreover, as it will subsequently be shown, t 
this low viscosity layer, as deduced from the phenomenological mechanical model, 
significantly increases the stresses in the crust produced Ly the primary tecton- 
ic and seismic motivating force and is instrumental toward elevating these stress- 
es to levels that can produce seismic events. 
(16,17) 

Quite recently, Tryggvason in an attempt to determine the characteristics of 

a model in which an elastic plate overlies a fluid substratum and which is consis- 


tent with the tilt measurements obtained in Icela:d, found that the viscosity of 


the fluid layer beneath the elastic plate cannot exceed 10!3 poise. This important 
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result obtained by Tryggvison based on tilt measurements from 1967 to 1971 fur- 


ther substantiates the basic nature of ou: phenomenolo:ical model of the earth, 
in so far as our fundamental phenomenological model inevitably led to enalytical 
investigations that, in fact, predicted the existence of such a low viscosity lay- 
(14,15) 
er in the upper mantle two years before Trygqvasun's results were publ ished. 
Another important aspect of this historical development concerns the concep- 
tion and mathematical formulation of a simple constitutive model of the crust and 
mantle which concurs with our basic phenomenological mechanical model and which 
quantifies, with approximation, the seismic process that led to the conception of 
this phenomenological model. This particular constitutive model provides a math- 
ematical framework that links analytically and numerically, the observable and 
monitorable features on the surface of the earth, such as displacements, strain, 
tilt, etc. measured in a region around an active fault with the attendant changes 
in the stress field within and throughout the crust and upper mantle, and with 
the relaxation time constants of the Maxwellian mantle. In this procedure, we 
adopt and follow a working principle based on a selective rule, namely, that in 
the construction of constitutive models designed to quantify the phenomenological 
description of the seismic mechanism, we first select the simplest constitutive é 
members of the class of Maxwellian materials, admitted by this description, and 
then progressively increase the complexity of these constitutive models only in 
so far as dictated by refined mzasurements of observable phenomena and by new facts. 
This, we believe, insures that the analysis «and numerical calculations are not 


overburdened by unnecessary complexity of the constit itive models for the purpose 


of accomodating the given data and furnishing the information sought. 


Within this matematical framework, a three-dimens ional boundary value prob- 


lem in static linear elasticity was formulated as an approximate formal charac- 

(18,19) 
terization of the phenomenological model for the crust-mantle system of the earth. 
Ths system is considered as a two-layered half-space, where elastic and inelastic 
Properties are attributed to both media. !t was assumed that both materials are 
incom:ressible. A Volterra dislocation, representing a fault plane, is imbedded 
in the upper layer, lying in a vertical plane with slip parallel to the free sur- 
face. The problem is solved by means of a modified form of the Galerkin vector 
and use of double Fourier transforms. Certain analytical solutions are obtained 
in quadrature for the surface displacements and for some stresses at the inter- 
face. Time dependent properties are introduced into the model by means of Biot's 
correspondence principle of linear viscoelasticity. The exhaustion of the in- 
formation rending capacity of this mathematical model of the seismic mechanism 
requires extensive numerical calculations which are readily tractable, but which 
we were unable to carry out because of lack of funding. 

We wish to stress, in this regard, that this mathematical model, if fully 
exploited, xipeeey provides a framework for developing a numerical model of the 
seismic mechanism which would numerically link the surface measurements, such 
as tilt, strain, etc., with the stresses within and throughout the crust and 
mantle. The development of such a numerical model of seismic mechanism has in- 
deed quite recently become a major concern of individuals actively committed to 


earthquake prediction. 
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11. THE CONCEPTUAL BACKGROUND OF THE PRIMARY TECTONIC AND SEISMIC MOTIVATING FORCE 


In the original phenomenological description of the seismic process which 
led to the conception of the phenomenological mechanical model of the Earth, as 
discussed in Part I, the source of the dynamical factor responsibie for motivat- 
ing the seismic process described therein was considered as time independent, but 
otherwise not specified nor restricted. However, a Primary Tectonic and Seismic 
Motivating Force (PTSMF) was subsecuently conceived by the senior author in 1965 
on the basis of the Universal Physical Constants, which led to the concept of Pri- 
macy of Force and the Principle of Maximum Uniformity. In this part of the report, 
we shall outline this previous development that provided the basis for the numeri- 


cal studies on verification of the PTSMF as subsequentiy car: ied out under this 


present contract. 


A. Theoretical Basis cf the PTSMF----Primacy of Force 


The focus and impetus which led to the conceptual framework in which the 
PTSMF was discovered are the Universal Physical Constants. This conceptual frame- 
work necessarily assigns to all forces in nature a position of primacy by ascrib- 
ing to them the role of constantly linking and everywhere dynamically reconciling 
the universal and immutable features of the time-space structure with specific 
phenomena and actual events, that is; the fundamental role of dynamically recon- 


ciling Universality and Specificity in all natural phenomena. 


* The discovery of this motivating force was outlined in an invited lectur«: 
presented at the National Earthquake Research Center of the United Stio3s 
Geological Survey, Menlo Park, California in March 1965, and subsequently 
amplified in lectures given at the Department of Geology, University of 
California at Davis and the National Earthquake Research Center of the 

United States Geological Survey, Menlo Park, California in 1973. 


This development necessarily compelled a critical examination of the Principies 


of Classical Mechanics with the expectation of being able to substantiate within 

its limited conceptual framework the inferences drawn on the primacy of force from 
the Universal Physical Constants. For this purpose, the senior author has reform- 
ulated the principles of classica! mechanics by adopting a procedure which treats 

force as primitive, that is; as fundamental and undefined. This leads to the con- 
cept of a force configuration which is then used to precisely define, on the basis 
of Newton's fourth axiom on the compounding ot forces, a particular class of force 
configurations, which are endowed with a particular dynamical aspect of uniformity, 
namely the class of equilibrium force configurations. This class of precisely de- 
fined equilibrium force configurations was then used to formulate what might appro- 
priately be called the fundamental principle of mechanics insofar as it evidently 

incorporates in a single proposition Galileo's principle as well as Newton's second 
and third laws. Moreover, its domain of validity is not necessarily restricted to 
classical mechanical phenomena insofar as this fundamental dynamical principle of 


force equilibrium is not predicted on the model which Galileo's principle and 


* 
Newton's law of motion are necessarily limited. This proposition asserts: ‘That 


every realizable force configuration is necessarity an equilibrium force conf ig- i” 


| uration and that this aspect of uniformity, namely eyuilibrium, is constantly, 


—— 


uniformly and precisely maintained in nature, at all times, at all locations and 
for all bodies.' From this formulation, it directly and incisively follows that: 
(a) The restrictions on mechanical phenomena expressed and/or implied by the 


fundamental principle of mechanics do not define and, in principle, do not deter- 


mine the forces in terms of which they are enunciated. Consequently, force is 


a SS 


See footnote on page 13. 
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not subsumed under the known principles of mechanics, since in principle, it is 
not determined by them. Therefore, the forces in terms of which the known prin- 
ciples of mechanics are formulated are escentially independent of these principles. 

As corollaries of (a) we obtain: 

(b) The principles of classical mechanics (PCM) do not determine the forces 
by which they are enunciated and by which they are dynamically and uniformly sus- 
tained in nature. Forces are therefore in principle not determined by the PCM; 

(c) That the PCM do not distinguish and consequently do not specify mechan- 
ical systems, either formally or actually; 

(d) That the PCM do not discriminate between conservative and non-conserv- 
ative mechenical systems; 

(e) That the PCM are aspecific. That is, the restrictions they uniformly 
invoke on all mechanical systems do not distinguish between them and consequently 
do not specify them, 

(f) That the PCM are histroically passive, that is; they neither deny or 
affirm the historical process, 

(g) That the PCM neither affirm nor deny reversibility; 

(h) That the PCM are acausal aid do not axiomatize the causal, i.e.; the 
evolutionary context of the forces by which they are enunciated. 

These developments produced the concept of mechanical speciation and the 
realization that forces, in addition to sustaining dynamically an aspect of uni- 

formity which, according to the fundamental law cited above is unifermly maintain- 
ed in all classes and subclasses of mechanical systems, also specify and thereby 
distinguish particular classes, as well as identify physically actua mechanical 


systems. Indeed, from these results we see that within the conceptual framework 


of classical mechanics, forces maintain the dynamical features which are universal 
for all mechanical systems, as well as specify particular classes and actual dynam- 
ical systems; thereby, in fact, linking and reconciling their universal and speci- 


fic mechanical features. 


Forces Unify and Specify. They Reconcile Unity in Diversity. 


B. The Discovery of the PTSMF 


Based on the Universal Physical Constants and the Principle of Maximum Uniform- 
ity, we conceived, in 1965, a primary tectonic and seismic motivating force (PTSMF). 
This primary tectonic and seismic motivating force is linked with the historical 
development of the earth and is the jynamical aspect of its mechanical evolution. 
However, it is not implied by the principles of classical mechanics since, as ex- 
plained in Section A, these principles are essentially historically passive, that 
is; they neither affirm nor deny the historical process on which this primary 
motivating force depends. All mechanical phenomena produced in the earth are accord- 
ingly understood to be essentially products of its mechanical evolution and the dy- 
namical consequences of the primary tectonic and seismic motivating force. 

The PTSMF essentially originates from six principal factors: 

(1) The rotation of the Earth; 

(2) The inelasticity of the Earth, as originally inferred from a phenomeno- 
logical description of the seismic process, a description which incorporates vir- 
tually indisputable phenomen logical features of thls process; 

(3) The Maxwellian upper mantle; 

(4) The non-uniformities in mechanical constitution of the Earth which spec- 


, ify and thus identify the actual Earth mechanical ly; 


(20-23) 
(5) The Principle of Maximum Uniformity applied conceptually and heuristic- 


ally to produce algorithms that link analytically and numerically the category of 
general information, which contains the fundamental principle of mechanics, to the 
categories of specific and concrete ian wate include answers to specific 
questions pertaining to the mechanical evolution of the Earth and which are not 
implied by the known propositions of classical mechanics; and 

(6) The development of specific hydrodynamical ramifications of these algor- 
ithms which are relevant to the calculation of stresses and deformation produced 
by the PTSMF in the mantle, and which have led to the formulation of the hydro- 
dynamical Variational Principles that augment the Navier-Stokes Equations. These 
include the original formulation of the Hydrodynamical Principle of Minimum Dissi- 
P| eee applied by Prigogene, et.al. to biophysical questions with some 
success, and the Hydrodynamical Principle of Virtual nigtncee used in conjuc- 
tion with the Finite Element soaee i a numerical augmentation to the Navier- 
Stokes Equations. 

Subsequent to the conception of the primary tectonic and seismic motivating 
force noted above, an analytical and numerical procedure was formulated for the 
purpose of calculating the stress and deformation fields produced by this force 
in mechanical models of the Earth. This analytical and numerical procedure was 
the basis of the original proposal that led to the initiation of this contract 
for the numerical verification of the PTSMF. Since the scientific and practical 
ramifications of this comprehensive development in Geomechanics remain to be 


fully explored, which unfortunately we have been unable to do under the present 


contract because of insufficient funding, the detailed presentation of this 


development will be given in subsequent publications when this task is satisfactor~ 


ily fulfilled. 


a 


111. NUMERICAL RESULTS AND VERIFICATION OF PTSMF 


A. Preliminary Remarks 


Observable consequences and seismic ramifications of the primary tectonic 
and seismic motivating force (PTSMF) are analytically and numerically examined by 
calculating the stresses and deformations it produces in a hierarchy of me-hanical 
models which constitutively portray and thereby specify the Earth with increasing 
mechanical non-uniformity. The results of these numerica' calculations thus pro- 
vide an important test of verification of the PTSMF and were produced with this 
principal objective in mind, namely, to ascertain whether the stress fields pro- 
duced by the PTSMF even in comparatively primitive mechanical models of the Earth 
can, in order of magnitude as well as in configuration, account for the stresses 
required to produce and to sustain tectonic and seismic activity in known seis~ 
mically active areas. For this purpose, it appears sufficient to consider limit- 
ing Maxwellian mechanical models of the Earth in which the relaxation times are 
virtually infinite. The stress and deformation fie'ds produced by the PTSMF in 
these limiting Maxwellian models may consequently be considered as snapshots of 
the stress and deformation fields evolving historically under this force in space 


and time. 
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a Our study of predictable and observable manifestations of the Primary Tectonic 


q and Seismic Motivating Force in terms of a hierarchy of mechanical models of the 

| Earth which progressively depict its mechanical non-uniformity with increasing spe- 
cificity, has led to a new basis for developing a deep and power ful synthesis of 

| the geological and geophysical perspectives on the science of the Earth. Particular 
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members of this hierarchy of mechanical models are distinguishable by the order 
and degree of specificity by which they depict the mechanical non-uniformity of 
the actual Earth. 

To fix ideas we find it helpful to introduce here the concept of orders of 
non-uniformity which pertain to fundamentally distinct aspects of the mechanical] 
non-uniformity of the Earth. 

By this concept, the prominent mechanical] non-uniformity of the Earth which 
pertains to i - globally stratified structure is conceived as a first order mech- 
inical non-uniformity which necessarily has global ramifications and which accounts 
for the production of a high threshold level of stress in the crust which origin- 
ates from the primary seismic motivating force. The high threshold of stress assoc- 
iated with this first order non-uniformity is an aspect of what might appropriately 
be called the primordial stress field of the Earth, since this field is associated 
with a mechanical non-uniformity which evidently was established during the very 
early history of the Earth and which has since persisted with considerable stabil- 
ity. The stability of this first order mechanical non-uniformity and of the assoc- 
iated primordial stress fie'd which it produces in conjunction with the Primary 
Tectonic and Seismic Motivating Force, is understood here to be an aspect of the 
Principle of Maximum Uniformity which is the theoretical basis of our work in geo- 
mechanics. 

Accordingly, the concentric layered material sphere is a configuration which 
reconciles within a finite domain maximum geometrical uniformity with minimum sub- 
Stantial non-uniformity. The homogeneous (unstratified) sphere is of course an 
ideal body which unifies within a finite domain the highest attainable aspects of 


substantial and geometrical uniformity, but which, however, is virtually unreal- 


izable when its parts are not uniformly subjected to motivating forces. When such 


mechanical non-uniformities prevail and are sustained by the physical environment 


of a homogeneous sphere, then according to the principle of maximum uniformity, 
the homogeneous sphere progressively adopts to the persisting mechanical non-uni- 
formities that are impressed upon it by producing a homogeneously concentric lay- 
ered structure in which the mechanical non-uniformities are localized at the inter- 
faces that join the virtually homogeneous iayers. By so doing, the concentric 
homogeneously layered sphere accomodates persisting mechanical non-uniformities 
which may be produced by its environment, by deviating minimumly from an ideally 
uniform, that is, homogeneous sphere. These observations help us to place in pro- 
per perspective the very encouraging results obtained and presented in this part 
of the report from a number of models belonging to the hierarchy of mechanical 
models of the farth, noted above. 

The first resu'ts were obtained for a set of mechanical models of the Earth 
which incorporate only the first order mechanical non-uniformities which obtain 
from the essentially concentrically stratified structure of the Earth, which ona 
global scale is highly stable. The primordial stress field produced by the Primary 
Tectonic and Seismic Motivating Force in conjunction with this first order and 
highly stable mechanical non-uniformity of the earth, ina significant sense, vir- 
tually dominates the stress field of the Earth, except where higher order mechanic- 
al non-uniformities which depend on latitude and longitude are comparatively signi- 
ficant. Such higher order mechanical non-uniformities are relatively significant 
in the crustal layer and particularly in transitions from oceanic to continental 
structure. 

In this context, we can refer to the mechanical non-uniformities associated 


with deviations in the primordial layered structure of the Earth from perfect 
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spherical symmetry, as mechanical non-uniformities of the second order. The mech- 


anical non-uniformities associated with changes in crustal structure and constitu- 
tion, that correspond to the oceanic and continental regions, we shall refer to as 
mechanical non-conformities of third order. The mechanical non-uniformities assoc- 
iated with increasingly local features which mechanically specify, and thereby 
indentify, a particular location of the Earth comprise a set of mechanical non- 
uniformities with increasingly higher order. 

We have calculated the stress fields produced by the Primary and Tectonic and 
Seismic Motivating Force in a set of models of the Earth which incorporate only 


, first order mechanical non-uniformities which are associated with the primordial 
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layered structure of the Earth depicted as a concentric homogeneously layered 


sphere. The results of this first set of calculations which, we believe, corres- 
pond to the primordial stress field of the Earth, are most encouraging, as they 
reveal prominent latitude dependent stress patterns which correspond to known 
tectonic features in areas which have been carefully explored. These areas in- 
b clude Anchorage, Alaska and California, among others. Moreover, the first set of 
calculations produced stress fields which in magnitude configuration and dis- 
tribution give compelling reason to ascribe the mechanical aspects of the evolu- 
tion of the Earth to the Primary Tectonic and Seismic Motivating Force, which was 
originally corceived in 1965 as the dynamical manifestation of its mechanical 
| evolution. The historical devlopment of the figure of the Earth and its tectonic 
and seismic features are in this context understood to be particular aspects of 


its mechanical evolution. 


r The results obtained in the second set of calculations were obtained for a set 


of models of the Earth which incorporate first order, as well as second order 
mechanical non-uniformities which include mechanical non-uniformities that corres- 
pond in first approximation to the oceanic and crustal features of the earth. 

These results imply tectonic features on a global scale which are strikingly sub- 
stantiated by observation. Specifically, these results imply a mirror image re- 
lation between the San Andreas Fault with its attendent stress field, and a corres- 
ponding prominent fault system at the latitude of the China mainland. Such a 

fault system has indeed been recently reported by Tuzo Wilson fellowing a recent 


visit to the Institute of Geophysics in China. 


B. Mechanical Models and Tests for Verification 
sneranean ee ae One UES ES Ol Verran On 


(1) Seto 
The first set of mechanical models in this sequence which we desig- 
nate as the zeroth set consists of homogeneous spheres which are consequently 
maximally uniform in geometry as well as :n mechanical constitution. The non- 
uniform force field created by the uniformrtation about a diameter of the sphere 
produces a non-uniform stress field with respect to latitude. However, the magni- 
tude of the stresses produced in this non-uniform stress field is insufficient to 
account, in order of magnitude, for stresses necessary to produce seismic phenomena 
Therefore, we shall not consider the results obtained in this case further. 
(2) Set t 
(a) Mechanical Non-uniformities 
The mechanical models of Set | incorporate the mechanical non- 
uniformities of primordal earth which were produced early in its history and evident- 
ly sustained since then with considerable stability. This set of mechanical mod- 
els consists of four concentric layered spheres representing the core, lower mantle, 
upper mantle and the crust. The mechanical properties, except the density, are 
assumed uniform within each layer, but differ from layer to layer. This set of ; 
models therefore preserves maximum uniformity in geometry, as in set zero, but has 
a minimal non-uniformity in its mechanical constitution. 

For a limiting Maxwellian model where the viscosity coefficients \' and yu! 
are virtually infinite, the mechanical properties can be characterized by the den- 
sity o and the two elastic moduli A and yp. Altogether, numerical calculations have 
been performed for eight different combinations of values of p, \, and pt assigned 


to each of the four layers with A'and u' infinite throughout (Table |) in this set 
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of mechanical models. The first seven models assume that the density varies in 
the radial a from the center of the core to the surface of the earth as 
shown in ne ae However, the density is assumed constant across each layer in 
model VIII (see Table ! and Figure 1). 
(b) The Numerical Results for the Stresses 

Models | to VI collectively demonstrate the effect of the shear 
modulus UW, as it varies from layer to layer, on the magnitude of the stresses. In 
general, it can be stated that as the value of y is progressively increased from 
the core to the crust in various models while \ and p are maintainec fixed in 
these models, higher values of the stress components a (in east-west direction) 


oo 


and Sp4 (in north-south direction) (Figure 2) are produced by the PTSMF in the 


upper layers. To illustrate this, let us consider the numerical result for O44 
obtained from models IV and VI (Figure 3). We note that the mechanical properties 
are the same in the core, in the lower mantle, and in the crust of these two models 
(see Table |), whereas the value of py in the upper mantle of model VI is about half 
the value assiqned to the upper mantle of model IV. As a result, the magnitude of 
Oo in the upper mantle of model VI is smaller than it is im model IV. We see, how- 


oo 


ever, that this greatly increases the stress O44 in the crust of model VI as com- * 
pared with model !V. This clearly shows that a reduction in the competence of a 
particular layer to support a Icad in shear produced by the PTSMF when it is ina 
fixed state of deformation results in an increase in level of stress in an overly- 
ing layer which has a higher shear modulus. Significantly, these particular numer- a 
ical results concur beautifully with the phenomenological description of the seis- 


mic process based on relaxation phenomena and which led to the conception of a 


phenome nological Maxwellian mechanical iodel of the Earth; a phenomenological model 


endowed with mechanical features that are essential to produce and sustain the 


PTSMF as an evolutionary mechanical force. 

it is further noted that stress concentrations in cen and Ogg are produced by 
the PTSMF at interfaces of two layers where discontinuity in the mechanical prop- 
erties occurs. These stress concentrations attain maxima at the interface joining 
the crust and upper mantle as well as at the interface joining the upper and low- 
er mantle. Moreover, in this set of primordal mechanical models, the normal stress 
component Co. (in the radial direction) is in general small. It is of the same 
order of magnitude as o,, and o in the core and lower mantle, but decreases rapid- 


oo 68 


ly toward the surface. in addition, the shear stress O. in these models is every- 


8 
where small. 
(c) Tests for the Verification of PTSMF 
Among the primordia] models 1! through VI, only models 1V through 
VI, and in particular model VI, yield sufficiently high stresses to cause tectonic 
and seismic events in the upper 12% of the Earth's surface; i.e., in the upper 
mantle and the crust. Consequently, only the numerical results obtained from mode] 
Vi will be used here to test the verification of the PTSMF as the primary motivator 
of the tectonic and seismic phenomena of the Earth, both globally and locally. 
The numerical results of the stress components C59 and Ou obtained from 
model VI are presented in Figures 4 and 5, respectively. As these results are ob- 
tained in the primordial model in which mechanical properties do not vary in the 


t-direction; i.e., with respect to longitude, the magnitude of these stress com- 


ponents remains unchanged for any particular latitude. O59 is, in general, com- 


pressive, except in the upper 2-3% of the crust at and near the equator. However, 


T9% changes gradually from compression in the polar region to tension at a lat- 
itude of about 45° and attains a maximum in tension at the equator. Since in this 
model De and Og are small everywhere in the mantle and the crust, eer and [99 
are essentially the principal stresses. The principal planes therefore are perpen- 
dicular to the sperical layers and are at right angles to the longitude and lati- 
tude. Consequently, the planes of maximum shear stress are at + 45° with respect 
to the longitude and latitude. 

In order to verify the hypothesis that the PTSMF is primarily responsible for 
motivating the tectonic and seismic phenomena of the Earth, both globally and local- 
ly, we shall examine the numerical results obtained for the stresses, in their mag- 
nitude and distribution, and resolve them along certain known geological features 
such as the San Andreas Fault in California, U. S. A., the Tancheng-Lukiang Fault 
Zone in Shantung, China and the region around Anchorage, Alaska, U.S.A. (Fiqures 
6 and 7). 

The San Andreas Fault is situated approximately between the latitudes 30° and 
4O°N with a strike in NW-SE direction. It is a recognized fact that the land mass 
west of the San Andreas Fault is moving in a northwesterly direction relative to 
the land mass east of the fault. This motion clearly implies a large shear stress 
along this fault, which is necessary to produce it and dynamically sustain it. 

The calculated primordial stress field in this region yields a tension in O66 and 
compression in Ig9° Since as noted above "46 and Op, are, for all practical pur- 
poses, the principal stresses, one of the planes of maximum shear stress is at 45° 
with the principal planes and therefore at that particular region coincides with 
the San Andreas Fault. Furthermore, the orientation of the shear stress along 


this plane coincides with the known movements of the land masses on the two sides 
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of the fault, i.e.; the shear stress is directed toward southeast on the eastern 


side and northwest on the western side of the fault; and the ma: nitude of this 
shear stress is sufficiently high and is equal to 165 bars; about fifty percent 
of the lower limit of the crushing strength of granite (see Figure 8). 
Correspondingly, the Tancheng-Luking Fault in ss dan runs in SW to NE 
direction from Shantunqg Province, Po Hai to Lianoning Province in the north (from 
30°N to SO°N latitude) may be considered, according to the description of the 
Chinese geophysicists as reported by Professor J. T. Wilson, as the mirror image 
of the San Andreas Fault produced with respect to the Pacific Ocean. In this re- 
aion, the land mass east of the fault is apparently moving in a northeasterly di- 
rection relative to the land mass west of the fault. Here again, one of the planes 
of maximum shear stress resulted from the calculated primordial stress field coin- 
cides with this fault and the orientation of the resolved shear stress is in accord 
with the movements of the land masses on the two sices of the fault. The magnitude 
of this shear stress is equal to 165 bars. 
Now, let us consider the region around Anchorage, Alaska, U. S. A.. Anchorage 
is situated at a latitude of about 60°N at the eastern end of the Aleutian Volcanic 
Arc. On March 27, 1964, an earthquake of magnitude 8.4 based upon surface wave ' 
amplitudes (M.) occurred in a vast region around Anchorage and major tectonic de- 
formation appeared in a segment of the continental margin lying between the Aleu- 
tian Trench and the Aleutian Volcanic Sa Post earthquake investigations suggest- 


ed that compressive stress tangential to the surface was dominant on land within ’ 


the region affected ity tectonic movements during the earthquake and the maximum 


-stress axis was oriented approximately normal, or at a small oblique angle to the 


major structural element within the focal region of the earthquake, ji.e.; the 


Aleutian Trench and the Aleutian Volcanic Arc. Resolution of the calculated pri- 
mordial stress field in this region yields a horizontal compressive stress equal 
to seventy-five percent of the breaking strength of granite in compression { 270 
bars) in the direction perpendicular to the Aleutian Volcanic Arc (Figure 8). 
This is in very good agreement with the results of the investigations mentioned 
above. 

In summary, the calculated primordial stress field obtained from certain 
mechanical models of the earth in Set | as a consequence of PTSMF has yielded a 
pattern of stress distribution which is in remarkable agreement with the known 
geological activities along the San Andreas Fault and the Tancheng-Lukiang Fault 
as well as in the region around the Aleutian Volcanic Arc and the Aleutian Trench. 
The magnitude of these calculated primordial stresses obtained from such mechan- 
ical models is sufficiently high for sustaining the tectonic and seismic activ- 
ities in such regions, if not for producing these features. From this, we con- 
clude that PTSMF is indeed primarily responsible for motivating the tectonic and 


seismic phenomena of the earth, both globally and locally. 


(d) Remarks .” 
Calculations have been performed to determine the effect of the i 
density and elastic modulus i on the calculated stress and deformation fields. In 
model VIII, the value of \ in the core has been increased to four times the value . 
of model IV; and in model VIII, the density is assumed to be constant in each lay- 1 


er. In both cases, the magnitude and the pettern of the calculated stresses and 


deformations are not significantly different than those obtained f:om model IV. 


However, these results must be considered as indicative rather than conclusive with 
respect to the effects of \ and 0 on the stress and deformation fields until more 
numerical results can be obtained for wider combinations of \ and 0 values. 
(3) Set UI 
(a) The Mechanical Non-Uniformitics 

In Set Il, the prominent mechanical non-uniformities on the sur- 
face of the Earth, i.e.; the oceans and the con:inents, are incorporated into the 
mechanical models considered. In order to simplify the analysis and numerical cal- 
culation, the mechanical properties p, \ and p in the crustal layer are assumed 


to vary in both r and 4 directions in the following manner: 


p=! + % (r) cos 24 
X= 1+ f, (r) cos 26 
v= 1+ £3 (F) cos 2 


These, therefore, depict models of equal area for the oceans and the continents 
without sharp discontinuities as in the actual earth, in topographic and in mechan- 
ical properties at the boundaries of the oceans and the continents, and may be ap- 
propriately called idealized oceanic and continental models. {t should be kept 
in mind that the results obtained from these models would only give qualitative 
indications of the effect of the non-uniformities of the mechanical Properties in 
¢ direction on the stress and deformation fields as compared with the actual Earth. 
(b) The Numerical Results for ihe Stresses 

In the numerical calculations, the mechanical properties in the 

core, lower mantle and upper mantle are assumed the same as those of model IV, 


Set 1. f(r), f2(r) and f3(r) are taken as second degree polynomials in r with 
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maximum ar minimum values occurring at the mid-point of the crustal layer. Three 
different sets of values have been assigned to f;(r), f(r) and f3(r) in the nu- 
merical calculations. However, the numerical results obtained for the stresses 
with [fi(r = 0.125 a i = . = -0. st i r ing. 

[ a. ane 5 and [fotr) di. Petry os 0.625 are most interesting 

The general pattern of the stress and deformation fields obtained in all 
three mechanical models considered is, in two respects, the same as that obtained 
in Set I: 

] The str c 
(1) tress components 0. , O14, Op,, and 
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the upper mantle and crust. Therefore, the tangential normal stress components, 


are, in general, small in 


Sas and Jogr are essentially the principal stresses there. 

(2) Pe is compressive in the polar region, becomes a tensile stress near 
a latitude approximately equal to 45°, and finally reaches a maximum value in 
tension at the equator for all longitude. 

However, although T59 is also compressive in the polar region as in the prim- 
ordial models of Set |, it might change into a tensile stress at some latitude or 
remain compressive throughout dependent on the longitude. For this reason, the 
resolution of the two principal stresses, 0,, and o 
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shear stress, sometimes accompanied by tension, extending to the equator where as 


in the planes of maximum 


in the primordial models of Set |, such a region is restricted to certain inter- 
mediate latitudes only. Furthermore, the magnitude of the stress components var- 
ies in the ¢ direction , as well as in other directions, due to the additional 

mechanical non-uniformities in the 6 direction included in this set of mechanical 


models. These two new factors, namely, the variation in the pattern of Cp9 and in 
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the magnitude of O44 and og in the ¢ direction, both contribute to produce quite 


clearly zones of alternately high or low magnitude stress fields. Consequently, 
seismically active and inactive regions in terms of the stress field can be differ- 
entiated in a global scale. 
({c) Tests for the Verification of PTSMF 

As in Set |, we will proceed, by using the numerical results 
of the stress fiela obtained in this set of mechanical models, to further verify 
the hypothesis that the PTSMF is indeed the primary motivator for the tectonic 
and seismic activity of the Earth, both globally and locally. For this purpose, 
¢ = 0 is assumed to locate at the mid-points of the Pacific Ocean. 

The resolution of the calculated stress field obtained from this set of 
mechanical models along the San Andreas Fault, the Tancheng-Lukiang Fault and 
around the Anchorage region yields similar results as those obtained from the 
corresponding models of Set 1, both in magnitude and pattern of the stresses 
(Figure 9). In addition, the calculated stress fields are given in Figure 9 
for the seismically active regions around Tokyo (pure shear stress), Kamchatka 
(large compression and shear) and Managua (large tension and shear). The mag- 
nitude of the calculated stresses in these regions are all sufficiently large 
to produce seismic events. Unfortunately, verification of the calculated stress 
fields in these latter regions is not possible due to lack of measurements of the ! 
actual stress fields. 

Figure 10 shows the resolution of the calculated stress field along the Mid- 
Atlantic Ridge and the West Coast of Norway, Britain and Portugal. Measurements 


(31) 


of the actual stress field in these regions have been made by Hast. The technique 


of Hast's measurement system was developed permitting the determination of 
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direction and absolute value of rock stress in bore holes directly. The results 
of the measurenents show that the stresses are generally horizontal (i.e.; tan- 
gential to the spherical surface) compressive and shear stresses along the west 
coast of Norway, Britain and Portugal and where the Mid-Atlantic Ridge meets with 
Iceland. Since the calculated stress field yields er and C9 as the principal 
Stresses near the surface of the Earth, therefore the calculated stresses act on 
the horizontal surface only. Furthermore, the resolved stresses along the west 
coast of Norway, Britian, Portugal and Iceland clearly show that the calculated 
stress field consists of compressive stresses accompanied by shear stresses ir 
these regions. These are in good agreement with Hast's measurement (Figure 10). 
Moreover, along the northern portion of the Mid-Atlantic Ridge near Iceland, the 
resolved calculated stress is compressive and accompanied by a shear stress; where- 
as along the southern portion of the ridge toward the equator, the rift system, 
the calculated stress is acting in the direction very nearly perpendicular to the 
rift with sufficiently large tensile stress necessary for producing the rift. 


These results, again are in accord with the observed geological features of these 


regions. 


IV. CONCLUDING REMARKS AND OUTLINE FOR FUTURE DEVELOPMENTS 


A. Concluding Remarks Concerning the PTSMF 

The numerical results presented in this report are obtained from models 
which portray the primordial mechanical non-uniformities of the Earth in terms of 
Maxwellian materials for the limiting cases of very long relaxation times. There- 
fore, the results obtained from these limiting models may be considered as snap- 
shots of the numerical results which would be obtained from the corresponding prim- 
ordial models of the Earth in which the stress and deformation fields are both space 
and time dependent in an evolutionary sense. 

Using these primordial models, we calculated the stress and deformation fields 
produced by the PTSMF operating on the primordial mechanical non-uniformities of 
the Earth depicted in these models. From these numerical calculations, we found 
that the PTSMF produces, for certian specific plausible mechanical non-uniformities 
of the Earth, stress fields of sufficient magnitude to activate and to sustain 
tectonic and seismic phenomena in the Earth, both globally and locally. The geo- 
metrical resolution of these calculated stresses along prominant geological features, 
such as the San Andreas Fault, Tancheng-Lukiang Fault, Aleutian Trench and Volcanic 
Arc and the Mid-Atlantic Ridge, shows that the calculated stresses, both in magni- 
tude and in configuration, can account fora large percentage of the stresses neces- 
sary to maintain or activate these features. In regions along the coast of Ice- 
land, Norway, Britain and Portugal where actual stress measurements have been made, 
the calculated stress fields are in essential agreement with the measurements. 


These results clearly provide substantial tests for the verification of the 
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hypothesis that the PTSMF is indeed the primary motivator of the tectonic and seis- 


mic phenomena of the Earth, both globally and locally. 


Furthermore, the results obtained from the mechanical models of Set I! show 
that the stresses would be localized according ‘9 the variation of the mechanical 
properties in the ® direction as well as in the r direction and thereby differ- 
entiating seismically active from seismically inactive regions. This convincing- 


ly demonstrates, we believe, that there exists a virtually unique relationship 


between the mechanical non-uniformities of the actual Earth and its siress field 


and that this unique relationship is indeed dynamically determined by the PTSMF. 

In this regard, we would like to mention here that a recent work by D. Pine 
(32,33) 

and J. Shaham concerned with estimating the elastic energy level of a reservoir 

sufficient to sustain seismic activity in conjunction with Chandler wobble and 

continental drift, is based on a semi-empirical phenomenologi~al theory which is 

specifically constructed to link the elastic energy of this reservoir with a few 

global parameters which may be directly determined from observation. In their 

theory, observational data are used to calculate certain basic coefficients which 

in turn are used to calculate an energy level of the reservoir. Then by apply- 

ing the elastic energy reduction" principle, which is in fact a particular aspect le 

of our principle of maximum uniformity, the authors indeed demonstrate that the 

polar drift must be westerly, as actually observed, in order to comply with this 

"reduction principle! used in conjunction with their semi-empirical theory for 


calculating the release of elastic energy from the reservoir. Accordingly, the i 


theory is thereby constructed, tested and supported on phenomenological grounds 


when this reduction principle, which is extraneous and thus not implied by known 
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principles of physical theory, is invoked. Although it satisfactorily reconciles | 
phenomenologically prescribed observational data with basic phenomenological as- 

pects of polar wandering, this does not as yet constitute a test of verification 

of their calculated energy level of the reservoir or account for its source. 

It is of interest to note that, although their semi-empirical theory does 
not provide a theoretical basis for calculating the total energy of the elastic 
energy reservoir with which it is principally concerned and which they, in iact 
calculated semi-empirically and thus phenomenologically from observational] data, 
it nevertheless tacitly includes certain essential features as well as a princi- 
ple which are, respectively, particular aspects of certain fundamental] features 
and a new fundamental principle which are necessary to produce and to sustain the 
PTSMF, which was previously indentified as the principal motivator of tectonic 
and seismic phenomena in the Earth. These include (1) geometrical and mechanical 
non-uniformities of the Earth, (2) inelasticity, (3) rotation, (4) phenomenolog- 
ical approach to geomechanics and (5) the principle of maximum uniformity. Fur- 
thermore, the energy of the ''elastic energy reservoir'' which, on the basis of 
our work, js evidently produced and sustained by the PTSMF, is shown by Pine 
Shaham, on rather convincing semi-empirical and phenomenological grounds, to be 
sufficiently high in energy content to account for most geophysical processes. 
This clearly amounts to a further test of verification cf the PTSMF. 

In summary, since the numerical results obtained so far lead to certain de- 4 
finitive facts and evidence which strongly support the conclusion that the PTSMF 


is indeed the primary motivator of tectonic and seismic phenomena in the Earth, 


and since this force places stringent and highly specific restriction on 


— = 
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the relationship between the actual mechanical constitution of the Earth and 
known geological features and geophysical facts, it is necessary to critically 
examine and thereby revise the conventional elastic seismological models of the 
Earth on the basis of the PTSMF and the fundamental phenomenological mechanical 
model of the Earth presented in Part |, which, in turn, points up the fundamental 


need for the developments of mathematical theories of inelastic seismology. 


B. Specific Objectives and Tasks for Future Development 


When the test of verification of the PTSMF has been substantially com- 
pleted by reconciling the seismic data used to obtain the conventional elastic 
seismological models, and if the results continue to be affirmative, clearly 
then it would convincingly follow that all known and observable geological fea- 
tures which are mechanical in origin must be necessarily attributable to this force 
and to its operation on specific mechanical non-uniformities of the Earth which 
together produce these features. Therefore, these known geological features, be 
they structural or morphological in character, provide, in this context, a rich 
source of information rendered by the PTSMF for obtaining the mechanical non- 
uniformities of the Earth which are specifically responsible for producing them. 
Furthermore, it is also clear that this method can be carried over to the objec- 
tive of finding new and various kinds of energy and mineral resources which would 
be correspondingly specified by their specific mechanical non-uniformities. 

Therefore, once this important step for obtaining scientific consensus is 
completed, the PTSMF can then be used with confidence as a fundamental source of 


information and procedure for retrieving new and additional informatlon on the 


mechanical non-uniformities of the Earth which mechanically Identify the Earth 
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with increasing specificity. For this purpose, an effort will be directed, in the 
next phase of our work, toward incorporating realistically in the mechanical models 
of the Earth, in as much detail as possible, all the known and observed geological 
features and mechanically non-uniform properties on the surface of the Earth, as 
well as in the crustal layer, upper and lower mantle, and the core. These geolog- 
ical features may include the overall land mass distributions, mountain ranges, 
sea fioor structures, oceanic ridges, rifts, faults and trenches with increasing 
detail, and the known non-uniform distribution of the mechanical and geometrical 
properties, such as thermal non-uniformities in the crust and upper mantle, as well 
as deviations from spherical uniformity, will be considered. In order to be able 
to retrieve highly specific information from such highly specific mechanical models 
of the Earth, sophisticated numerical methods which have a'ready been developed 
and applied and which are based on a hydrodynamical variation principle and the 
finite element fon ae be used. This analytical and numerical procedure, aug~- 
mented by concurrent studies pertaining to the development of new methods in ex- 
perimental seismology which are based on new mathematical developments in inelas- 
tic seismology, will produce an economically feasible capability for retrieving 
very specific information on the mechanical non-uniformities of the Earth with 
highly practical and scientific ramifications. 

In order to effectively pursue the numerous and far reaching implications 
and scientific and practical ramifications of the PTSMF, noted above, in con- 
junction with the phenomenological mechanical model described in Part | of this 
report, we propose to initiate specific investigations pertaining to two principal 
objectives, namely, (1) earthquake prediction, control and modification and (2) 
the identification, with high specificity, of the mechanical non-uniformities of 


the Earth and its applications to energy and mineral exploration. These two 


ho 


principal objectives respectively concern (a) a systematic development of a hier- 
archy of mathematical theories and associated numerical methods of inelastic seis- 
mology and (b) the analytical and numerical investigations of a class of constit- 
utive models of the phenomenologically described seismic process that originally 
led to the conception of the basic phenomenological model; a model which has con- 
sistently anticipated and circumscribed subsequent developments in Geomechanics. 

A particular constitutive model of ‘he seismic process was originally describ- 
ed and mathematically treated tee Its formulation in present form already 
provides a basis for rationally developing a numerical model of the seismic pro- 
cess which operates in accordance with our original phenomenological description 
of the seismic process. This formulation links the surface displacements with the 
stress and deformation fields within the crust and mantle in regions surrounding 
the faults. Furthermore, it facilitates the numerical representation of progres- 
sively refined descriptions of this process as dictated by new phenomenological 
facts and thereby accomodates the representation of increasingly complex Maxwel |~- 
ian rheological features as demanded by these refined descriptions. Moreover, in 
this regard, we are now able to model and treat numerically the interaction of a 
crustal layer with materials of the upper mantle which displays non-linear hydro- ’ 
dynamical features arising from the convective transport of momentum in the under- 

(27) 
lying mantle. 
A fundamental aspect of the seismic mechanism as constitutively modeled in 


our earlier work is of course its cyclicity which depends upon the interaction of 


this constitutive model with a prevailing ambient stress field which was posited 


there, but without mathematical and numerical representation. With the discovery 


=. 


hy 


of the PTSMF and its subsequent verification, it is now possible to deal with this 
problem of cyclicity mathematically and numerically by bringing this force to bear 
on the seismic mechanism as represented by the hierarchy of constitutive anu mathe- 
matical models of the seismic process. This development, we believe, will also 
lead, as a by-product, to a rationale for understanding and interpreting observ-~- 
(34-36) 
able precursory indices, such as the V/s, anomaly, which may be useful for the 
anticipation of forthcoming earthquakes. 

As a further adjunct of this development, the capability concerning ear th- 
quake prediction and control can produce information that will be useful for ob- 
taining new discriminating indices for identifying artificial earthquakes that are 
mechanically designed to camouflage its detection. This necessarily requires re- 
fined physical understanding of the production of wave trains and of their trans- 
mission in non-uniform inelastic media and in their associated amtient stress 
fields, since such knowledge is required by an adversary to mechanically camou- 
flage the event. Such new discriminating indices must therfore ultimately be 
traced to the fine differences of the phenomenological mechanical features that 
attend the production of the natural and artificial earthquakes and are consequent-~ 
ly responsible for the subtle, but yet discernable, differences ir the so-called 
quality of the seismic signals emanating from them. This is analogous to the unambig- 


uous differences in the quality of production well known in musical acoustics. 


C. Further Observations Concerning Inelastic Seismology 


This section refers to the marked and radical change produced in theo- 


retical seismology with the advent of the realization <hat the inetasticity of the 


Earth is a fundamental and thus an essential aspect cf its seismicity. 
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Theoretical sesimology started essentially as a mathematical theory of elastic 
sesimology with a basic committment to retrieve from measurements certain informa- 
tion about the mechanical constitution of the Earth; namely, the distribution of 
its elastic and inertial properties. With this limited objective in mind, Jeffreys, 
among others, persuasively demonstrated that the retrieval of this particular mech- 
anical information about the Earth, from seismic data and from constitutive mech- 
anical models, was insensitive to the modeling of the inelasticity of the Earth. 
This position was developed from the standpoint of a comparatively limited per- 
spective on the seismicity of the Earth, which evidently did not consider the 
most fundamental mechanical aspect of its seismicity, namely the mechanical phe- 
nomena that attend the evolution and the recurrence of earthquakes. 

Thus, theoretical seismology evidently started with the preconceived idea 
that the elasticity and the inertia of the Earth are its only fundamental mech- 
anical properties and accordingly devised mathematical seismological models which 
are SOlely committed to the identlflcation of these two properties withIn the Earth. 

But, if one starts out as we did with the realization that inelasticity is a 


mechanical property on which the seismicity of the Earth depends fundamental ly, 


then clearly the commitment to its numerical indentification in terms of appro- tes 

\. 
priate inelastic constitutive models of the Earth used in conjunction with seis- y 
mic data is virtually inescapable. It is this important new development which 


is essentially based on the comparatively recent discovery that the inelasticity 
of the Earth is fundamental to its seismicity, that marked and initiated the rad- 
ical change in theoretical seismology; a change which leads to the emergence of 


mathematical inelastic seismology. 
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However, the transition from elastic to inelastic mathematical seismology | 
presents a formidable challenge in the selection of appropriate constitutive in- 
elastic models from an indefinitely large number of classes of conceivable in- 
elastic materials, a challenge that does not present itself in elastic seismology 
which refers strictly to only one class of materials. To grapple with this dif- 
ficult task, and as a matter of procedure, we refer again to the phenomenological 
mechanical model of the Earth inferred from the phenomenological description of 
the seismic process which dominates the class of all admissible inelastic con- 
stitutive models. By doing so, the indefinitely large number of classes of con- 
ceivable inelastic materials is reduced to one admissible class, namely the class 
of Maxwellian materials. Here again we see that the phenomenological mechanical 
model indeed provides definite and stringent guidelines and criteria for selec- 


ting a class of admissible inelastic constitutive models. 


D. Concluding Remark 


As a concluding remark, we note that the historical description given 
here of our work in geomechanics, extending over a period of more than two de- 
cades, shows that this work is consistently substantiated by results subsequently 
obtained by other investigators. This important fact clearly demonstrates the 


power of and constantly adds credibility to our basic capability as a fundamental 


source of information in geomechanics. 
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Figure 2. Configuration of the Stress Field in Spherical Coordinate System 
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Figure 3. Distributions cf the Calculated Primordial Stress Field 4 
with Depth for Different Values of 2 
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Figure 4. Distributions of the Calculated Primordial Stress Field g6 with Depth 
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Figure 5. Distributions of the Calculated Primordial Stress Field O44 with Depth 
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Appendix A 


Numerical Results Obtained for 


Model IV, Set | 
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